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ABSTRACT 

We report observations of type III radio bursts at decimeter wavelengths (type Illdm bursts) - sig- 
natures of suprathermal electron beams propagating in the low corona - using the new technique of 
radio dynamic imaging spectroscopy provided by the recently upgraded Karl G. Jansky Very Large 
Array (VLA). For the first time, type Illdm bursts were imaged with high time and frequency reso- 
lution over a broad frequency band, allowing electron beam trajectories in the corona to be deduced. 
Together with simultaneous hard X-ray (HXR) and extreme ultraviolet (EUV) observations, we show 
these beams emanate from an energy release site located in the low corona at a height below ~ 15 Mm, 
and propagate along a bundle of discrete magnetic loops upward into the corona. Our observations 
enable direct measurements of the plasma density along the magnetic loops, and allow us to constrain 
the diameter of these loops to be less than 100 km. These over-dense and ultra-thin loops reveal the 
fundamentally fibrous structure of the Sun's corona. The impulsive nature of the electron beams, their 
accessibility to different magnetic field lines, and the detailed structure of the magnetic release site 
revealed by the radio observations indicate that the localized energy release is highly fragmentary in 
time and space, supporting a bursty reconnection model that involves secondary magnetic structures 
for magnetic energy release and particle acceleration. 

Subject headings: Sun: flares - Sun: radio radiation - Sun: corona - Sun: magnetic topology 



1. INTRODUCTION 

Certain types of solar activity - flares and jets - are 
powered by the impulsive release of energy through fast 
magnetic reconnection. In the standard scenario, a re- 
conn e ction region located i n low corona ( Masuda et al.l 
119941 : iSui fc Holmanl 120031 : iBattaglia fe Be nz 2006|) ac- 
celerates particles to nonthermal energies. These par- 
ticles carry a large fraction of the released energy and 
play an important r ole in energy transport processes 
(jLin fc Hudson! fl 9 76[ ). Electron beams can form through 
the propagation process, if not from the acceleration 
mechanism itself (jBastian et al.l Il998). and travel up- 
ward and downward in the corona along magnetic field 



lines. The downward propagating beams will eventually 
be stopped by the dense chromospheric plasma and pro- 
duce hard X-ray (HXR) emission. The chromospheric 
plasma is heated by these beams to high temperatures 
and expands hydrodynamically into the corona where it 
i s observed i n EUV and/or soft X-ray (SXR) wavelengths 
(Benz 2008). However, direct information about these 
fast electron beams and the region from whence they 
originate has remained elusive. 

One method to probe these beams is through radio 
observations. The radio signature of coronal electron 
beams, first discovered in the lat e- 1940s at meter wave- 
lengths (jWild fc McCreadvlll950h . is called a type III ra- 
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dio burst. It is the result of the nonlinear conversion of 
Langmuir waves generated by beam instabilities to elec- 
tromagnetic radiation at the fundamental or harmonic 
of the local electron plasma frequency, f pe « 9^Jn~ e kHz, 
where n e is the electron number density. Since n e varies 
with height in the corona, the radio emission from the 
electron beam drifts from high to low frequencies for 
upward-propagating beams, and in the opposite sense 
for downward-propagating beams. Their high frequency 
counterpart - type Illdm bursts - w as first observed 
in the early 1960s (jYoung et al.l ll961). and was subse- 
quently found to be cl osely associated with impulsive 
energy release in flares jKane 1981; Denn is et al.l 11984] : 
lAschwanden et all 119851 IWhite et al.l 1201 lh . For these 
reasons, type Illdm bursts are believed to be an impor- 
tant diagnostic of impulsive magnetic energy release. 

A barrier to exploiting type Illdm emissions has 
been the lack of true "radio dynamic imaging spec- 
troscopy" , where simultaneous imaging observations are 
available at each frequency and time that a burst is ob- 
served. Imaging observations of type III, type Illdm, 
or type U bursts (a type III in a closed magnetic loop) 
have been reported at one, o r a fe w, discrete frequen- 
cies (e.g.. lAschwanden et al.l 119921: iRaulin et al.l 119961 : 
lAurass fc Kleinlll997l : iPaesold et al.ll2001h . but radio dy- 
namic imaging spectroscopy has not been possible. Here, 
we report the first use of radio dynamic imaging spec- 
troscopy with the upgraded VLA to observe type Illdm 
bursts during a coronal jet event that accompanied a so- 
lar flare. In §2 we discuss the observational setup of the 
VLA and the data. We also present observations from 
a number of complementary instruments. In §3 we dis- 
cuss the implications of these data and place them in an 
interpretive context. We briefly conclude in §4. 

2. OBSERVATIONS 

First dedicated in 1980, the VLA has been recently 
upgraded with state-of-art receivers an d electronics 
which have greatly increased its capabilities (jPerlev et al.l 
l201lh . While the VLA is a general-purpose radio tele- 
scope, provisions were made to enable solar observing. 
Particularly important are the relatively large instan- 
taneous bandwidths for imaging, the large number of 
spectral channels, and the fast sampling times available, 
which enable the new observing technique of radio dy- 
namic imaging spectroscopy. 

The first solar observations by the upgraded VLA were 
made on 2011 November 5 when the VLA was in the 
D configuration, for which the longest antenna baseline 
is 1 km. The observations were made in the 1-2 GHz 
frequency band (A = 15 — 30 cm) using 1024 spectral 
channels, each 1 MHz in bandwidth, with a time resolu- 
tion of 100 ms. Due to limitations in data throughput 
at that time, only 17 antennas could be used. These 
were nevertheless sufficient for imaging the source of ra- 
dio emission. The angular resolution of the 17- antenna 
array was 35" x 95" at 1.5 GHz at the time of observa- 
tion, which scales linearly with wavelength. Both senses 
of circular polarization were observed. Unfortunately, a 
hardware error resulted in the loss of most of the data 
from 1.5-2.0 GHz, rendering it unsuitable for imaging al- 
though dynamic spectroscopic information is available. 
We therefore focus on the data from 1.0-1.5 GHz. Cal- 
ibration of the antenna gains and frequency bandpasses 



was performed by referencing the observations to stan- 
dard sidereal sources. When observing the Sun, 20 dB 
attenuators were inserted in the signal path. Their pres- 
ence introduces fixed perturbations to the antenna gains 
that were measured prior to the observations and subse- 
quently corrected by calibration. 

Type Illdm bursts were observed in association with a 
coronal jet during the aftermath of a GOES class Ml. 8 
SXR flare (Fig. IB) in active region (AR) 11339. The 
cross-power dynamic spectra are shown in Fig. 1G-H, 
where each bright vertical or nearly- vertical feature rep- 
resents an individual burst. The largest group of bursts 
occurred between 21:21:40-21:22:10 UT (Fig. 1H), well 
correlated with a peak in the 12-25 keV HXR emis- 
sion observed by the Ramat y High Energy S olar Spec- 
troscopic Imager (RHESSI) (|Lin et al.l 120021 ) (Fig. II, 
the red curve), indicating the release of nonthermal elec- 
trons. At 21:21:00 UT, another (weaker) HXR peak can 
be seen, during which the radio continuum has a small 
enhancement. A small cluster of narrow-band type Ill- 
like bursts can be found around the time of this HXR 
peak (Fig. 1G). The coronal jet was observed by the So- 
lar Dynamics Obs ervatory (SDO) At mospheric Imaging 
Assembly (AIA) (|Lemen et al.ll2012l ) (Fig. 1A) during 
this energy release event, followed by subsequent bright- 
ening of closed loops (Fig. 1C-F). The jet was also ob- 
ser ved by the EU V Imager (EUVI) aboard STEREO 
B (|Wuelser et al.ll2004f ) that observed the Sun from the 
east (Fig. 2). Interestingly, no concurrent interplanetary 
type III bursts were d etected by the STEREO/WAVES 
(|Bougeret et al.ll2008[ ). 

Dynamic imaging spectroscopy with the VLA allows 
each pixel in the dynamic spectrum to be imaged. Al- 
though the type Illdm bursts themselves were unresolved 
by the VLA, the radio source centroid position can be 
determined with high accuracy. The high signal-to-noise 
ratio of these bursts allows their positions to be fit to 
~ 1% of the nominal resolution, or < 1 Mm. Exam- 
ples of the fitted centroid locations for two temporally- 
resolved bursts (indicated by the arrows "1" and "2" 
in Fig. 1H) are shown in Fig. 3B-C and Fig. 3D- 
G, colored from blue to red corresponding to increas- 
ing frequency. At any given time, the centroid locations 
at different frequencies follow well-defined trajectories, 
distributed from southwest (bottom-right) to northeast 
(upper-left) with increasing frequencies. All the trajecto- 
ries fall within an envelope with the high frequency end 
(high density, thus low coronal heights) originating near 
the location of the EUV jet and the HXR sources low in 
the atmosphere (Fig. 3A, showing a superposition of all 
the type-IIIdm centroid locations from 21:20:30-21:22:10 
UT). In addition, for bursts that are temporally-resolved, 
the emission drifts from high to low frequencies with a 
rate of 0.3-1 GHz/s, and their centroid locations move 
to the southwest direction with time, indicating that the 
electron beams were propagating upward away from the 
EUV and HXR source (e.g., Fig. 3D-G). These observa- 
tions show the electron beam trajectories directly and, 
hence, the coronal magnetic field lines along which they 
streamed away from the reconnect ion site. 

3. RESULTS AND DISCUSSION 

The close spatiotemporal association of the type- 
Illdm-burst trajectories and the HXR footpoints (Fig. 
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Fig. 1. (A) Composite image of AR 11339 showing emission in 304 A (red, 0.1 MK), 94 A (green, 6 MK), and 171 A (blue, 0.8 MK) 
from SDO/AlA. (B) GOES SXR flux showing the Ml. 8 flare (peaks around 20:35 UT) and its aftermath. The coronal jet occurred between 
21:18-21:25 UT (dashed vertical lines). (C-F) Dynamic evolution of the coronal jet and the subsequent loop brighting in AIA 131 A (0.4 
MK and 10 MK) images. The field of view (FOV) is indicated by the large box in Fig. 1A. (G-H) VLA dynamic spectra in which type 
Illdm bursts (the bright vertical or nearly-vertical features) are present. (I) VLA radio (blue) and RHESSI 12-25 keV HXR (red) light 
curves from 21:20:30-21:22:30 UT, showing the temporal correlation of the radio flux and type Illdm bursts with the HXR flux. The small 
and large boxes correspond to the time ranges of the dynamic spectra shown in (G) and (H). 
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Fig. 2. — The coronal jet (marked by the white arrows) seen 
by both STEREO/EUVI B (left column) and SDO/AIA (right 
column) at around 21:26 UT. SDO was viewing AR 11339 from 
the Earth, while the STEREO B satellite was trailing the Earth 
by 103° at the time of the observation thereby viewing AR 11339 
from the east. X- and Y-axes are aligned with solar east- west 
and south-north respectively in each perspective. White crosses 
mark the same location at the foot of the EUV jet, based on direct 
coordinate transformations. 



3B-G) suggests that the X-ray-producing downward- 
propagating electron beams and the type-IIIdm-emitting 
upward-propagating electron beams originate from 
a common energy release site, confirming previous 
ideas based largely on the temp oral correlation be- 
tween HXR/type-IIIdm emission (|Dennis et al.l Il984t 
lAschwanden et all Il995h . It also places the reconnec- 
tion site at a location above the HXR footpoints and 
below the height of the highest-frequency type-IIIdm- 
burst sources; the latter is estimated to be ~ 15 Mm as 
we discuss further below. This energy release event was 
likely associated with evolving satellite magnetic polari- 
ties in AR 11339 (Fig. 4A-D). They coincided spatially 
with the HXR footpoint source (Fig. 4A-D) and the 
base of the EUV jet (Fig. 4E-H). The evolution of the 
HXR source and the EUV jet suggests that energy re- 
lease and plasma heating took place in two stages: first, a 
HXR footpoint source appeared over the northern satel- 
lite polarity (Fig. 4A) during the first HXR enhance- 
ment around 21:21:00 UT, while an EUV jet was initi- 
ated (Fig. 4F; the timing was estimated within the 12-s 
cadence of the SDO/AIA images) that involved eruption 
of an arcade-shaped feature. During the second HXR en- 
hancement peaked around 21:21:50 UT, the HXR foot- 
point source moved gradually to the southern satellite 
polarity (Fig. 4B-D), while another eruptive EUV jet 
was initiated (Fig. 4G-H). These multi- wavelength ob- 
servations allow us to propose a self-consistent scenario 
describing the physical processes that occurred during 
this event (Fig. 5). First, interaction of the satellite 
magnetic polarities with the surrounding magnetic flux 
triggered magnetic reconnection above the northern neu- 
tral line around the first HXR peak (indicated by the "X" 
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Fig. 3. — (A) Emission centroids of all type Illdm bursts observed from 21:20:30-21:22:10 UT, colored from blue to red in increasing 
frequencies, showing electron beam trajectories in projection. Background is the SDO/AIA 131 A image at 21:22:09 UT. Red contours 
are the 12-s integrated 12-25 keV HXR emission during the second HXR peak (around 21:21:50 UT). Emission centroids of a temporally 
resolved type Illdm burst observed from 21:21:52.4-21:21:52.6 UT (shown by arrow "1" in Fig. 1H) are shown by (B-C) for two successive 
100-ms integrations. Another temporally-resolved type Illdm burst from 21:21:57.1-21:21:57.5 UT (arrow "2" in Fig. 1H) is shown by (D-G) 
for four successive 100-ms integrations. Red contours are the 12-25 keV HXR source with a 4-s integration closest to the type-IIIdm-burst 
times. 



symbol) . Both upward- and downward-propagating elec- 
tron beams were produced, resulting in the observed type 
Illdm bursts and HXR footpoint emission (the shaded 
area). The released magnetic energy can trigger instabil- 
ities and initiate "blowout " EUV jets that in volve erup- 
tion of a magnetic arcade ( M oore et al.ll2010l ). Magnetic 
field lines were stretched and led to further reconnec- 
tions above the southern magnetic neutral line as well. 
The accelerated particles led to subsequent brightenings 
of the southern EUV loops, HXR footpoints and type 
Illdm bursts. 

The VLA data can be used to derive the electron num- 
ber density n e along the type-IIIdm-emitting loops since 
it is directly related to the plasma frequency or its har- 
monic. We assume that the emission is likely harmonic 
plasma radiation for two reasons. First, fundamental 
plasma radiation tends to be more highly circu lar ly- 
polarized than its harmonic (|Dulk fc Suzukil [1980) . Yet 
the majority of the observed type Illdm bursts is weakly 
polarized (<20%). Second, fundamental plasma radia- 
tion results in a group delay of > 0.5 s in its transit to 
Earth compared to HXR photons, while th e delay is ex- 
pecte d to be far less for harmonic radiation (|Dennis et al.l 
11984( 1. No such delay > 0.5 s is found when the HXR 
and radio data are cross-correlated. Therefore, we have 
n e « (//18 kHz) 2 = 3.3-7x 10 9 cm -3 over the frequency 
of 1.0-1.5 GHz, distributed over a height range of 30-80 
Mm for the majority of the type Illdm source centroids. 
The radio source height was obtained by measuring the 
projected distance from the HXR footpoints along the 
beam trajectories. Based on an estimate of the inclina- 
tion angle of the EUV jet (~45°) from the simultaneous 



SDO/AIA and STEREO/EUVI B stereoscopic observa- 
tions (Fig. 2), and assuming a similar inclination for the 
type Illdm bursts, the de-projected height is obtained. 
Fitting the ensemble of type Illdm trajectories to a sim- 
ple hydrostatic density model yields a best-fit density 
scale height L n = n e {—dn e ldK)~ x « 40 Mm, correspond- 
ing to a temperature of 0.8 MK assuming hydrostatic 
equilibrium. Returning to the location of the reconnec- 
tion site, the height of the highest-frequency type-IIIdm- 
burst sources is found to be as low as ~ 20 Mm from the 
imaging results. Considering that some bursts occurred 
at frequencies as high as 2 GHz, for which imaging is 
not available, the maximum height of the reconnect ion 
site can be reduced to ~ 15 Mm within the context of 
this density model. A striking feature of the type-IIIdm- 
source distributions is their frequency variation across 
the envelope of trajectories (Fig. 3A); that is, for a 
given frequency the emission is not necessarily located 
at the same height from burst to burst, suggesting that 
the type-IIIdm-emitting loops are inhomogeneous per- 
pendicular to the magnetic field, with An e /(n e ) « 25%. 
For the limited number of temporally-resolved bursts the 
electron beam velocity can be directly measured from 
the movement of the type Illdm source centroids in time 
(e.g., Fig. 3D-G), which gives Vb ~ 0.3c, where c is the 
speed of light. This is consistent with the electron energy 
responsible for the observed HXR emission in the energy 
range 12-25 keV, further demonstrating that their parent 
electrons are accelerated from a common site. 

Surprisingly, no trace of loop-like structures can be 
found along the type-IIIdm-burst trajectories against 
the background in any SDO/AIA EUV filters (cf. 
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Fig. 4. — (A-D): A series of 12-25 keV HXR emission (red contours) overlaid on the SDO Helioseismic and Magnetic Imager (HMI) 
(|Scherrer et al.l \201% magnetogram observed during the two HXR peaks around 21:21:00 UT (A) and 21:21:50 UT (B-D). (E-H): A 
detailed view of the coronal jet showing its temporal evolution in SDO/AIA 131 A images. The FOV of the images is indicated by the 
small box in Fig. 1A. 




Fig. 5. — Schematic illustration of the magnetic field evolution, 
reconnection, EUV jet, and generation of the electron beams and 
HXR footpoint sources during the impulsive energy release event. 
See the text for descriptions. 



iRaulin et al.lll996[ ), implying that the column emission 
measure £ of the type-IIIdm-emitting loops is too small 
to result in detectable emission or absorption relative 
to the background. To quantify this, we evaluate the 
observed intensity of an AIA filter image I(x, y) = 
£(Tb, x, y)R(To)At, where £(Tq, x, y) is the emission mea- 
sure assuming an isothermal plasma temperature To, 
R(To) is the filter response function at To and At is 
the exposure time. With the estimated temperature (0.8 
MK) of the type-IIIdm-emitting loops, the most proba- 
ble EUV band to detect any signature is the AIA 171 
A band. A negative result means the intensity contri- 
bution AI(x,y) by these loops is smaller than the rms 
intensity variation 07. In regions where the type Illdm 
source is present we measure 07 ~ 40 — 100 DN pixel -1 
for AIA 171 A images with At = 2 s. With R 171 « 10" 24 



DN cm 5 s _1 pixel -1 at 0.8 MK (jLemen et al.ll2012l ). the 
upper limit of the emission measure of these loops is es- 
timated to be £ max = Gjj AtjRyjx < 5 x 10 25 cm -5 . 
With £ = n 2 V/w 2 = n 2 Nd 2 /w, where V is the to- 
tal volume of N identical loops along the line of sight 
within the instrument resolution element of size w ( 0.4 
Mm for AIA), the upper limit of the loop diameter is 

tax = (U^/^e) 1/2 /^ 1/2 = l/N 1/2 . For a density of 
n e = 3.3 — 7 x 10 9 cm -3 inside the loops, I is found to 
be only a few xlO to ^100 km under rather generous as- 
sumptions, and d max would be even smaller with a filling 
factor 1/N < 1. 

The background plasma density and temperature can 
be estimated by utilizing the observed background in- 
tensities Ib(x,y) from all the AIA EUV bands sensi- 
tive to a wide ran g e of coronal temperature. We use 
lAschwanden et all (|2011[ Ts automated forward- fitting 
technique to obtain a differential emission measure distri- 
bution, d^(T, x, y)/dT, which yields an integrated emis- 
sion measure of £5 ~ 5 x 10 2T cm -5 peaking at ~2.8 MK. 
With £5 = f n e (h) 2 dh, the background density is found 
to vary from 6.9 x 10 8 to 4.8 x 10 8 cm -3 over 30-80 Mm 
at which the type Illdm bursts occur. Hence the type- 
IIIdm-emitting loops may be cooler (by a factor of ~3) 
and denser (by roughly an order of magnitude) relative to 
the background coronal m edium. This i s con sistent with 
the scenario proposed by iBenz etaLl (|1992f ) for allow- 
ing the type Illdm emission to escape. The over-density 
is also consistent with past conclusions based on the 
observed contrast between £ toward bright EUV/SXR 
loops and 6, toward the background corona (|Hara et al.l 
1 19921 : lAschwanden! l2QQlh , which suggests that the rele- 
vant active region loops may have already been supplied 
with high-density chromospheric plasmas through previ- 
ous heating processes. 
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The existence of ultra-fine magnetic structures in the 
coron a has been previously suspected (see, e.g. iReald 
|2Q1Q[ and references therein) but no current instrumen- 
tation can resolve such structures at coronal tempera- 
tures. Recent observations of ultra-fine magnetic struc- 
tures h ave been of plasma at chromospheric temper- 
ature s (|Ji et al.l 120121 ; lAntolin fc Rouppe van der Voortl 
|2012[ ), reporting loop diameters ranging from ^ 100 to 
a few xlOO km. Our observations suggest the corona 
is "fibrous" in nature, consisting of many unresolved 
"strands" . It should also be noted that the type-IIIdm- 
emitting strands have a direct connectivity to the re- 
connection region. Thus these strands should be even 
thinner down at the reconnection region with a diameter 
of at most 10s of km due to their general expansion with 
increasing coronal height (DeForest 2007). Furthermore, 
the multitudes of discrete electrons beams observed have 
access to spatially distinct coronal magnetic strands in 
<1 s (Fig. 3A), which indicates that the reconnection 
region likely consists of a large number of discrete re- 
connection sites in a localized spatial volume. Our ob- 
servations rather directly suggest a bursty reconnection 
scenario involving a localized reconnection region con- 
taining a distribution of many small-scale dynamically- 
evolving structures for magnetic energy release and par- 
ticle acceleration. This picture is qualitatively consis- 
tent with previous wor k on fragmentary e nergy release 
on the Sun (see, e.g. Ivan den Oordl 11994 and refer- 
ences therein), as well as more recent developments in 
the physics of magnetic reconnection, such a s experi- 
ments in the labora tory ([M oser & Bellan 2012) and nu- 
merical simulations (jKliem et al. 2000; Dr ake et al .1120061 : 
iKarlickv fc Ba rta 200?]), which show the reconnection re- 
gion consists of many "magnetic islands" or "filaments" 
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developed through magnetohydrodynamic instabilities 
such as tearing or kink instability. The constraints in 
the temporal and spatial scales of fragmentations we re- 
port here and in future observations of this kind may 
help guide theoretical reconnection models as well as con- 
strain parametric inputs to numerical simulations. 

4. CONCLUSIONS 

Dynamic imaging spectroscopy of type Illdm bursts 
with the VLA has allowed us to map the trajectories 
of electron beams produced by magnetic energy release 
during a coronal jet that occurred in the aftermath of a 
flare. Electrons escaped along discrete, ultra-fine strands 
into the upper atmosphere, producing the observed type 
Illdm bursts. Together with magnetic, EUV, and HXR 
data, we showed that the energy release site was associ- 
ated with a coronal jet and that the magnetic reconnec- 
tion process involved different locations at different times 
and that, furthermore, the density at each location var- 
ied significantly. The spatial scales in the reconnection 
sites are likely 10s of km or less. We conclude that the 
magnetic energy release process is highly fragmentary 
and that the surrounding coronal medium is fibrous in 
nature. 
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